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Spin- Triplet Superconductivity Mediated by Phonons 
in Quasi-One-Dimensional Systems 

Yuuichi SUGINISHI and Hiroshi Shimahara 

Department of Quantum Matter Science, ADSM, Hiroshima University, 
Higashi-Hiroshima 739-8530, Japan 

We investigate the spin-triplet superconductivity mediated by phonons in quasi-one- 
dimensional (QID) systems with open Fermi surfaces. We obtain the ground state phase di- 
agrams. It is found that spin-triplet superconductivity occurs for weak screening and strong 
on-site Coulomb interaction, even in the absence of any additional nonphonon pairing inter- 
actions. We find that the nodeless spin-triplet state is more favorable than the spin-triplet 
state with line nodes, for the parameter values of the QID superconductors (TMTSF)2X. We 
also find that QID open Fermi surface, which is the specific feature of this system, plays an 
essential role in the pairing symmetry. We discuss the compatibility of the present results with 
the experimental results in these compounds. 

KEYWORDS: anisotropic superconductivity, spin-triplet pairing, phonon-mediated pairing interactions, 
quasi-one-dimensional systems, organic superconductors, (TMTSF)2X 



1. Introduction 

The relation between the symmetry of the supercon- 
ducting order parameter and the mechanism of the pair- 
ing interaction in exotic superconductors is one of the 
most important subjects in condensed matter physics. 
In this paper, we are interested in this subject in quasi- 
one-dimensional (QID) systems, in connection with the 
organic superconductors (TMTSF)2X with X = PFe, 
CIO4, etc. 

The compound (TMTSF)2PF6 exhibits superconduc- 
tivity in the proximity of spin-density-wave (SDW) in the 
phase diagram.^ Hence, at an early stage after the dis- 
covery, the possibiUty of the superconductivity induced 
by antiferromagnetic spin fluctuations was proposed.^ 
Takigawa et al. have measured the NMR relaxation rate 
Tf ^ in superconducting (TMTSF)2C104.^ Their results 
suggest the anisotropic superconducting order parame- 
ter with line nodes on the Fermi surface as Hasegawa 
and Fukuyama argued.^ These experimental and theo- 
retical results on T^^ seem to be consistent with the 
mechanism of the superconductivity induced by the an- 
tiferromagnetic spin fluctuations. In this mechanism, the 
phase diagram of SDW and d-wave superconductivity, 
whose order parameter has line nodes, has been repro- 
duced semiquantitatively.^^^ 

However, recent experimental results seem contra- 
dictory to one another with respect to the symmetry 
of the order parameter. Bclin and Behnia obtained a 
temperature dependence of the thermal conductivity in 
(TMTSF)2C104 that indicated a nodeless superconduct- 
ing order parameter.'' Their result seems inconsistent 
with the NMR result mentioned above, which indicates 
the line nodes. As a possible consistent explanation of 
these results of the NMR and the thermal conductiv- 
ity, a nodeless d-wave superconductivity mediated by 
antiferromagnetic spin fluctuations has been proposed.^ 
Furthermore, in (TMTSF)2PF6 at P ~ 6 kbar and 
(TMTSF)2C104, the upper critical field Hc2iT) seems 



to exceed the Pauli paramagnetic limit Hp at low tem- 



perature 



9, 10 



From this result of Hc2, several authors 



have discussed possibilities of spin-triplet superconduc- 



tivity. 



However, there has not been any observation 



of a reentrant superconducting phase, which should oc- 
cur in spin-triplet superconductors at higher field. 
The experimental result Hc2 > ^^p can be explained by 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO or LOFF) 
state^^' and strong coupling effects. Recently, Lee et 
al. have measured a Knight shift in (TMTSF)2PF6 at 
P ~ 0.65 GPa, the absence of which strongly suggests 
spin-triplet pairing . ^ ^ 

Here, we shall briefly summarize the experimental re- 
sults and pairing symmetries. The full-gap state without 
any alternation of the sign, which we call s-wave state, is 
consistent with the thermal conductivity, but could not 
explain the absence of the Hebel-Slichter peak of T-^^ at 
all,'^ as well as the power low behavior and the Knight 
shift. It is also diflicult to explain the phase diagram. 
Therefore, the s-wave state is most unlikely, at least at 
the present. The singlet state with four line nodes, which 
we call the d-wave state, appeared the most likely at an 
early stage from the phase diagram and the NMR re- 
sults, but it could not explain the recent Knight shift re- 
sults, unless any high-field triplet state is assumed. The 
thermal conductivity may be explained even in the d- 
wave state with anion order.® The spin-triplet states are 
compatible with the absence of the Knight shift. For the 
open Fermi surface, both the full gap spin-triplet state 
without any line node and that with two line nodes at 
ky — are possible. We call the former and latter states 
the px-wave and pj,-wave states, respectively. However, 
the data of the NMR T^^ and the thermal conductivity 
are diflicult to be explained consistently. The p^^-wave 
state could explain the thermal conductivity, but could 
not explain the temperature dependence of T^^ cx 
for Tc/2 <T <Tc. The coherence peak becomes smaller 
than that of the s-wave state,'' but there has not been 
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any satisfactory theoretical fitting at the present. The 
Py-wave state could explain the line node, but it is dif- 
ficult to reproduce the thermal conductivity. In the Px 
and Py-wave states, the pressure dependence of Tc does 
not appear to be easily explained as in the d-wave state. 
Therefore, there is not any perfect explanation consistent 
for all experimental results at the present. 

Many theoretical studies have also been presented in 
order to clarify the mechanism of superconductivity in 
these compounds. In QID Hubbard models, consistently 
with the early result in the random phase approximation 
(RPA),'' the fluctuation exchange (FLEX) approxima- 
tion by Kino and Kontani^* and the third order pertur- 
bation theory by Nomura and Yamada^^ have showed 
that d-wave superconductivity is favorable. Spin-triplet 
/-wave-like pairing was also examined by Kuroki et al?^ 
and by Tanaka and Kuroki. There might be a consis- 
tent explanation on the basis of /-wave pairing, although 
/-wave pairing interactions arc usually weak, because 
the /-wave states oscillate many times in the momen- 
tum space. Recently, Kohmoto and Sato have discussed 
the possibility of nodeless spin-triplet superconductivity 
in a QID system with electron-phonon interactions and 
antiferromagnetic spin fluctuations.^^ 

The anisotropic pairing does not necessarily mean 
the nonphonon mechanism of the pairing interactions. 
Anisotropic pairing interactions mediated by phonons 
have been proposed by several authors. Foulkes and Gy- 
orffy argued that p-wave pairing is favorable in metals 
such as Rh, W and Pd, if the short range Coulomb in- 
teraction suppresses the s-wave pairing. In connection 
with high-Tc cuprates, Abrikosov obtained an anisotropic 
s-wave order parameter in a model of the phonon- 
mediated interaction with weak Coulomb screening .24.25 
Moreover, Friedel and Kohmoto, and Chang et al. have 
shown that rf-wave superconductivity is induced in a sim- 
ilar model with additional interactions mediated by an- 
tiferromagnetic spin fluctuations. 2^ Varelogiannis has 
also studied phonon-mediated d-wave superconductivity 
in K-BEDT-TTF compounds. More recently, one of the 
authors and Kohmoto proposed that anisotropic super- 
conductivity can be induced by the phonon-mediated in- 
teraction between electrons in ferromagnetic compounds, 
in connection with UGe2.^^ They also found that the 
anisotropic pairing interactions are enhanced in layered 
compounds with a large layer spacing.'^" In their pa- 
per, they discussed that spin-triplet superconductivity in 
the compound Sr2Ru04 and the organic superconductors 
may be induced by the phonon-mediated interaction. 

In the compound (TMTSF)2C104, very large isotope 
shift of Tc was observed, which suggests the large con- 
tribution of the phonon-mediated interaction to the su- 
perconductivity. The large isotope shift may be explained 
by considering the spin fluctuations.^^ 

In this paper, we investigate the pairing symmetry of 
QID systems within the phonon-mediated pairing inter- 
actions. Recently, one of the authors have derived effec- 
tive Hamiltonians of anisotropic superconductors from a 
model of electron-phonon systems with coexisting short- 
range and long-range Coulomb interactions.^^ We ap- 
ply these models to the QID system. In particular, we 



examine the phonon-mediated pairing interactions with 
screened electron-phonon interactions. Since we are moti- 
vated by the above controversy in the QID organic com- 
pounds (TMTSF)2X, we are particularly interested in 
which pairing symmetry is favored in the phonon mecha- 
nism for the parameter values of these compounds. How- 
ever, the range of our study is not limited to these c;oni- 
pounds. We are mainly interested in the effects of the 
open Fermi surface, which is one of the speciflc features 
of the QID systems. We show that nodeless spin-triplet 
state (px-wave state) is more favorable than the p^-wave 
state and the d-wave state in the phonon mechanism. 

This paper is constructed as follows: In the next 
section we introduce the formulation and examine the 
anisotropic superconductivity mediated by phonons in 
QID systems. In § 3, we show the results of the numeri- 
cal calculation and ground state phase diagrams. The last 
section devoted to the summary and brief discussion. 

2. Formalism 

2.1 Effective pairing interactions 

First, we consider a conventional form of pairing inter- 
actions mediated by phonons 



Ves{q,uJ + i5) 



-901 



(1) 



\q\^+qliv^ - [w(q)]2' 

where and w(q) denote the range of the pairing in- 
teraction and the renormalized phonon dispersion energy, 
respectively. We do not consider the detail of the mecha- 
nism of the electron-ion interaction, but only assume the 
range of the pairing interaction . In this sense, our 
model is phenomenological. This simplification is justi- 
fied when we are interested in a behavior of the pairing 
interaction in a length scale larger than the lattice con- 
stant. For shorter length scale comparable to the lattice 
constant, the lattice anisotropy also contributes to the 
anisotropy of the electron-phonon interactions. However, 
we neglect it for simplicity, since we are mainly interested 
in the specific feature of the QID open Fermi surface, 
which is taken into account through the momentum de- 
pendence in VeS by setting q = k — k' , where k and k' 
denote momenta on the QID Fermi surface. As is referred 
in ref. 34, near the Fermi surface where \uj\ < Lu{q), the 
interaction becomes attractive due to the overscreening 
effect. Except this, we ignore the dynamical effects in this 
paper. Within the weak coupling theory, the frequency 
dependence of the pairing interaction Ves{q,to + i6) is 
simplified as 



Ves{q) = -g- 



q^ + q. 



2 ' 



(2) 



for ICfcM^fc'l < ^D, and Vcs{q) = otherwise, where 
wd and denote the Debye frequency and the electron 
energy measured from the Fermi energy, respectively. We 
assume QID band structure with the dispersion relation 



= -2ta cos(A;^as) - 2tb cos{kybs) - n, 



(3) 



where /i, Og and feg denote the chemical potential, the 
intermolecular distances along the a and h directions, 
respectively. In eq. (2), we have introduced the effec- 
tive coupling constant g, which depends on the phonon 
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dispersion. Since it is difficult to estimate go and g 
from the first principle, we regard g as a parameter in 
this paper. Equation (2) has been studied by many au- 

thorS.24-27,29,30 

In the application of eq. (2) to the lattice system, it 
is more precise to replace \q\ with minjf \q — K\.^^ The 
effective interaction of eq. (2) has the peak around q = 0. 
Actually, the same structures should exist near q = K, 
where K denotes the reciprocal lattice vector. Hence, we 
obtain the effective Hamiltonian 



- q max - — 
K \q 



It 



(4) 



We refer to this Hamiltonian as model (a) in this paper. 
Secondly, we introduce a model to examine the efi'ect 

of the corrections due to the charge fluctuations. Within 
the RPA,^^''^^ the effective interaction is written as 



Ves{q) = -ffmax- 



■"K^\q-K\-^ + q^,Xo{q)IN{Qy ^' 

We refer to the Hamiltonian of eq. (5) as model (b). Here, 
Xo{q) denotes the static free susceptibility 



Xo{q) = 



(6) 



where f{x) = l/(c^/'^=^^ + 1). In the original form of this 
model, we must set qs = qTF- However, we extend the 
model by replacing q^p with a correct screening length 
q~^ so that it includes model (a) with qg ^ ^tf as a 
limiting case Xo('Z)/-^(0) ~ 1. This replacement of gxF 
by qs corresponds to taking into account the vertex cor- 
rection T{k,k',q) of Vq by an average quantity f with 

Thirdly, we introduce a model which includes the cor- 
rections due to the short-range Coulomb interaction U 
other than the direct repulsion. Since the details of the 
derivation arc presented in the previous paper, '^^ we shall 
describe only the outline here. We recall that the pairing 
interaction mediated by plionons should be proportional 
to M^/Ky{q) in the absence of U, where Mq and Ky{q) 
denote the electron-ion coupling constant and the dielec- 
tric function due to the long-range Coulomb interaction 
Vq, respectively. The electron-ion interactions are, more 
or less, the interactions between charges, irrespectively of 
their details. Therefore, in the presence of U, it is plausi- 
ble to replace Mq with Mq/nu{q), where Ku{q) denotes 
the dielectric function due to U. Furthermore, the po- 
larization function Xo{q) in the dielectric function 
should be replaced by that with a correction due to U, 
which we write xu{q)- In the RPA with respect to U, we 
have Ku{q) = 1 -|- Uxoiq) and 



. _ Xo(g) xo(g) 
l + UxQ{q) Ku{q) 



(7) 



Therefore, eq. (5) is modified as 



T4ff(g) 



9 



■ max 



[Mq)? ^ \q-K\^ + qMQ)IN(Oy' 



(8) 

which we call model (c) in this paper. 

The range of the effective interaction q~^ is on the or- 



der of the screening length in eq. (4), because it should re- 
flect the range of the Coulomb interaction between elec- 
trons and ions. In fact, if we adopt the Thomas- Fermi ap- 
proximation, we obtain the effective Hamiltonian of the 
same form as eq. (4) with q^ = ^tf- However, we should 
be careful when we use the Thomas-Fermi approximation 
for quantitative purpose on the screening length. For ex- 
ample, in the (TMTSF)2X systems, the Thomas-Fermi 
screening length is estimated as — 1.05 A, smce 
ta ~ 0.25 eV and yccii ~ 357.2 A^, where V^eii denotes the 
unit cell volume. ^^'^^ Hence, q^^ becomes much smaller 
than the intermolecular distance = 3.649 A in the a di- 
rection. However, it is obviously underestimation due to 
the Thomas-Fermi approximation. In practice, we could 
not apply the Thomas-Fermi approximation to phenom- 
ena of such a short length scale. In the real materials, 
the screening charges could not exist so densely within 
the distance shorter than the intermolecular distance ttg. 
Therefore, for the physical reason, it is more appropriate 
to assume q~^ > as rather than to put 5g~^ = 5xf ^^e 
application to the (TMTSF)2X systems. 

2.2 Gap equation 
The gap equation of superconductivity is expressed as 



(9) 



where we have defined Ek' = y £.k' + l^(^')P> and N 
denotes the number of the lattice sites. In this gap equa- 
tion (9) , the summation over k' in the right-hand side is 
restricted to | < lou near the Fermi surface. In the nu- 
merical calculations, we neglect the temperature depen- 
dence on Xo{q) since k^T <C ta - Further, we consider only 
intra-layer pairing and omit the A;2-dependence in A(fe) 
in this paper. Therefore, kz integral of Ves{q) appears 
in the gap equation. However, we have confirmed that 
the kz integral of Ves{q) has a peak around qx = qy = 
which can be fit in practice by a Lorenzian function with 
qs adjusted. Hence, we neglect g^-dependence in Vf.s{q) 
for simplicity. 

Now, we consider the quasi-one-dimensional (QID) 
systems. Once we have obtained the gap equation, the 
inter-plane electron hopping energy tc is negligible, while 
it justifies the BCS mean field approximation for suffi- 
ciently low temperatures. Near the Fermi surface, we put 
k ~ {±kFx{ky), ky), where the signs -|- and — correspond 
to the areas of the Fermi surface of fc^ > and kx < 0, 
respectively. Here, k-px{ky) is the Fermi wave vector of 
the kx direction at ky. Then, we can write 



V{k - k') ~ V{kFx{ky) T kFxik'y)^ - Kj) 

^\V^ + + HkyMy) = V^-){ky,k!y) 

W(+-)(fc„A:') = y(-+)(fc„A:') 



(10) 



and A(fc) ~ /:\'^^^^'^'^^''^\ky) . In the weak-coupling limit, 
one readily finds that the gap equation (9) near the su- 
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perconducting transition temperature takes the form 
1 r'^ bdk'„ 

w/b 



M(+)(fc,A_ _i r/' 



27r 



fV( + + \ky,k'^) V( + -\ky,k'y)\ /A( + )(fc;) 

\V(- + \ky,k'y) V^-\ky,k'y)) \A(-\k'y)j ' 



(11) 



where we have defined the density of states p{ky) 
l/4:iTtaSm[kFx{ky)a] and 



A 



= ln 



nksTc 



(12) 



with the Eiiler constant 7 = 0.57721 • • • . 

This matrix equation (11) is easily diagonaHzed as 

1 r^'' bdk 



with 



1 r'^/'> hc\k' 

A^(fc,) = -\j V^iky, k'y) A^ik'y) (13) 

{ky)-^P[ky) [^M^ky)-A(-){ky), 



A^iky) 

A%ky) 



(14) 



and 

V [Ky,Ky)-\^ ^ ^(++) _ , 

(15) 

where we have defined 

V(^''\ky,k'y) ^ ,f{ky)V^^^'Hky,k'y),/^). 

We expand the effective interaction and the order pa- 
rameter on the Fermi surface as 



00 CXJ 

V^'''\ky,k'y) = J2_12l^^nhm{ky)jn{k'y) 

+ V^n^lm{ky)Uk'y) 



m=0 n=0 



(17) 



and 



(18) 



m=0 



with a = s or a, which correspond to symmetric and an- 
tisymmetric states with respect to kx , respectively. Here, 
we have defined Jm{ky) = nmCOs{mkyb) and "ym{ky) = 
rim sm{mkyb) , and the normahzation factors rim = \/2 
for m 7^ and rim = 1 for m = 0. The expansion factors 
are given by 



Tr/b 
-7r/b 



bdky r/*" bdk'y 



27r 



TT/b 



2tt 



XJmiky)V^^^'Hky,k'y)^n{k'y) 



'I'' bdky r/"- bdk'y 



(19) 



n/b 



27r 



TT/b 



2tt 



X-frniky)V^^^'Hky,k'y)%{k'y). 

By inserting eqs. (17) and (18) into the gap equation 



Table I. The classification of the possible pairing states and rel- 
evant components of the gap functions. 





7m (fey) 


7m (fey) 


a = s 


Aq : «-wave state 


A| : py-wave state 




A| : d-wave state 




a = a 


Aq : pa;-wave state 


Aj : dxy-wave state 



(13), we obtain matrix equations 
00 

EA" A" = AA° 



(20) 



E 

n=0 



\a AO! _ \ AO 



where 



.('y(++) _ 

V mn mn / ' 

.(v(++) + v'-+-A, 

\ mn ' mn / ' 

V mn * mn / 



(21) 



We summarize our notation in Table 1. 

With the maximum eigenvalue A, the superconducting 
transition temperature Tc is given by 



2e'^ 

ksTc = wd exp [-1/A] , 



(22) 



for A > 0. 

Here, we note that the Coulomb pseudo-potential /Uq 
defined by 

UN{0) 

'^^ ~ 1 + UN{0) ln{Wc/ujB) ^^^^ 

must be subtracted from the eigenvalue A for the s-wave 
state obtained above. Here, Wc is the cutoff energy 
on the order of the bandwidth. It is obvious that the 
on-site Coulomb repulsive interaction U is not included 
in models (a)-(c), except the corrections due to U in the 
electron-phonon interaction and the screening function in 
model (c). However, it reduces only the s-wave pairing 
interaction, because of the symmetry. For example, in 
(TMTSF)2X, it is estimated as [/ ~ 1.5ia from the SDW 
transition temperature at tb = O.lta-^ 

At this point, wc have three parameters, the effec- 
tive electron-phonon coupling constant g, the screening 
length and the on-site Coulomb interaction U. Fur- 
thermore, we shall consider the electron number n per a 
site and a spin, the transfer integral tj, and the ratio ag/bg 
as given parameters within the QID systems with open 
Fermi surfaces. Here, it is easily verified by exchanging 
particle and hole that the results for n are equivalent to 
those for n — 1. For example, n = 1/4 and n = 3/4 give 
the same results below. 

3. Results 

In this section, wc show the results of the numeri- 
cal calculations. In Figs. 1-9, we set lvd = 200 K,^'^ 
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Fig. 1. Dimensionless coupling constants A as functions of the screening length ■ The labels (a), (b) and (c) correspond to the 
names of the models. The solid, dotted, dot-dashed and 2-dot-dashed curves show the coupling constants A for Px, Py, d and djij^-wave 
pairing, respectively. In the inset, the dashed curve shows the result for the s-wave state. 




Fig. 2. Dimensionless coupling constants A as functions of the ratio as/be. The labels (a), (b) and (c) correspond to the names of the 
models. The solid, dotted, dashed, dot-dashed and 2-dot-dashed curves represent the results of A for Px, Py, s, d and dxy-wave pairing, 
respectively. The Coulomb parameter is subtracted from A of s-wave pairing. 



tb/ta = 0.1 (except for Fig. 6), n = 1/4 (except for 
Fig. 5) and a^/hs = 0.468 (except for Figs. 2 and 7), 
considering the compounds (TMTSF)2X. These com- 
pounds have a slight strain of the lattice, but we have 
confirmed by the numerical calculations that it changes 
the result only very slightly. Therefore, we only show 
the results when the crystal axes are orthogonal. We set 
Wc/ton = 25.12 in /x^ from the experimental values,^'' 
since Wc = [W - fi') fi' ^ VSta for 1/4-fillcd band, 
where W and /z' denote the band width and the chem- 
ical potential measured from the bottom of the band, 
respectively. 

In Fig. 1, the results of the dimensionless coupling con- 
stants A are shown as functions of the screening length 
q~^. The labels (a)-(c) of figures correspond to the name 
of the models, which have been defined in the previous 
section. The s-wave coupling constant decreases rapidly 
with increase of the screening length (qsas)~^ and is close 
to anisotropic coupling constants in {qsas)~^ > 1. In a 
whole region of {qsas)~^, the Px-wa.ve state is more favor- 
able than the py-wscvc state. Especially for large {qsas)~^, 
the Px-wave state can easily overcome the s-wave state 
with an assistance of fi^. We find the similar behavior in 
models (a)-(c) for (^sas)"^ > 1. 



Figure 2 shows the dimensionless coupling constants A 
as functions of the ratio Og/bg. We assume (ggOs)"^ = 1, 
U/ta = 1.5 and g/ta = 5 as an example. The py-wavc 
state is favored for ag/bs > 1, while the p^-wave state is 
favored for Cs/bs < 1, when s-wave state is suppressed. 
In Figs. 1 and 2, it is also found by comparing the results 
of models (a)-(c) that the s-wave pairing interaction is 
weakened by the corrections due to the charge fluctuation 
and the short-range Coulomb interaction in addition to 

In Figs. 3-9, we concentrate ourselves to model (c). The 
results are qualitatively the same as those in models (a) 
and (b). 

Figure 3 shows the phase diagram at T = in the 

{qsas)~^ - U/ta plane. The parameter g/ta = 5 is taken 
as an example. It is found that the Px-vfave supercon- 
ductivity occurs in the region where the screening effect 
is weak and the short-range interaction U is strong. 

Figure 4 shows the phase diagrams at T = in the g/ta 
-U/ta plane. In each phase diagram, as is printed in it, 
the parameter [qsas)~^ = 0.3, 1 are taken as examples. 
It is found that the p^-wave superconductivity is favored 
for weak electron-phonon coupling and weak Coulomb 
screening. 




Fig. 3. Phase diagram in the ^-U plane at T = 0. Here, SC 
stands for superconductivity. 



Figure 5 shows the phase diagram at T = in the 
n - U/ta plane. We show the results in the region of n 
where the Fermi surface is open. We set {qsas)~^ — 1 and 
fl/ia = 5 as an example. The p^-wave superconductivity 
is favored for large n. 

Figure 6 shows the phase diagram at T = in the 
tb/ta - U/ta plane. We set (^sOs)"^ = 1 and g/ta = 5 
as an example. It is found that when th/ta increases, the 
Pa^-wave superconductivity is suppressed. Figure 7 shows 
the phase diagrams at T = in the Qs/bs - U/ta plane. 
The parameters, g/ta = 5, (qsOs)"^ = 0.3 or 1 are taken 
as examples. It is found that the Pa;-wave state is favored 
for small Og, while the Py-wave state is favored for large 
as- 

These results are physically interpreted as follows. The 
phonon-mediated pairing interactions are attractive in 
the whole momentum space of q = k — k', and have 
a peak at q = 0. Thus, the order parameter A(fe) at 
k is favorable if A(A;')'s of fe' « fe have the same sign 
as A(fc), while unfavorable if A(fe')'s at fe' w fc have 
the opposite sign as A(A;). Therefore, the py-wa.vc state 
is originally more unfavorable than the px-wiwc state, 
because of the node at ky = 0, where A(fc) changes its 
sign. However, when the areas of the Fermi surfaces of 
kx > and kx < approach by increasing Ug, the px-wave 
state becomes unfavorable, because A.{kx > 0,ky) and 
A.[kx < O.ky) have the opposite signs. The dependences 
on tb of the coupling constants are also explained in the 
same manner as the dependences on the ratio as /6s- Near 
the ends of the first Brillouin zone, the distance between 
the sheets of the Fermi surface becomes narrower as tb 
increase. (For tb/ta > 0.35, the Fermi surface is closed.) 
For the px-wavc state, the interaction becomes repulsive 
for pair hopping A(A;) A(fe') near the zone ends. As a 
result, as tb increases, s-wave superconductivity becomes 
more favorable than px-'wavc one. 

Figures 8 and 9 show the calculated momentum depen- 
dences of the order parameters. We assume (gsOs)"^ = 1 
and U/ta = 1.5 as an example in both figures. It is con- 
firmed that the py and dxywaye states have a line node 



Wq,a,= l.O 



p^-v/scwe SC 




Fig. 4. Phase diagrams in the g-U plane at T = 0. 



on the Fermi surface at fc^ = 0, the d-wave state has line 
nodes at ky w ±7r/26, and the s and p^-wave states are 
full gap states. 

4. Summary and Discussion 

Now, we summarize the results. We have examined 
pairing interactions mediated by phonons in quasi-one- 
dimensional (QID) systems with open Fermi surfaces. 
It was found that spin-triplet superconductivity occurs 
for weak Coulomb screening and strong on-site Coulomb 
interaction U, even in the absence of any additional non- 
phonon pairing interactions. We have examined two pos- 
sible spin-triplet states. One has a nodeless gap function 
A(fc) for the open Fermi surface, while another has line 
nodes a,t ky = on the Fermi surface. We call the former 
Pa;-wave state, while the latter pywave state. We have 
found that when the Fermi surface is open, as the ratio 
as/&s increase, p^-wave superconductivity becomes more 
favorable. We have explained the physical reason for it. 
It was argued that the Px-waYC state occurs when the 
conductive chains are separated with a sufiiciently large 
spacing. 

In (TMTSF)2X, the ratio a^/bs = 0.468 and the hole 
density n = 1/4 give sufficiently large separation of the 
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Fig. 5. Phase diagram in the n-U plane at T = 0. 
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Fig. 7. Phase diagrams in the ae/ba-U plane at T = 0. 



Fig. 6. Phase diagram in the tb-U plane at T = 0. 



fcj: > and fca; < sheets of the open Fermi surface 
when tb = O.lta- Hence, the p^-wave state is more fa- 
vorable than the py-wave state in those compounds. It 
was also found that the p^-wave state is more favorable 
than the 5-wave state for n = 1/4 and U ~ IMa, which 
are appropriate to (TMTSF)2X, if we set {qsas)~^ > 1, 
which is physically plausible. 

In our model, tb may be regarded as an effective pa- 
rameter that reflects an effect of pressure. If ag/bs and tf, 
increase with increase of the pressure, the coupling con- 
stant for p^-wave pairing decreases as found in Figs. 2 
and 6. This pressure dependence may contribute to the 
observed pressure dependence of to some extent. 

In order to discuss the reality of the phonon-mediated 
spin-triplet superconductivity, wc crudely estimate the 
parameters for (TMTSF)2X from the observed transition 
temperature ~ 1 K. Here, we assume p^-wave pairing. 
If wc assume cjd — 200 K"^^ and Tc ~ 1 K, we have 
A ~ 0.184. Then we find g ~ AMta for (gsOs)-^ ~ 1 
from Fig. 1. For such a choice of parameter values, in 
order to suppress 5-wave pairing, the on-site Coulomb 



repulsion must heU > 2ta — 0.5W, which seems realistic 
as the order of the magnitude. This value is rather larger 
than the estimation U « l.Sta from the SDW transition. 
However, it can be explained by taking into account that 
the repulsive on-site interaction is enhanced by the spin 
fluctuation near the SDW phase from the bare value U 
in practice. 

In conclusion, the nodeless spin-triplet superconduc- 
tivity is favorable in the compounds (TMTSF)2X, if s- 
wave superconductivity is suppressed. The present result 
consistent with the experimental data of the Knight shift 
and thermal conductivity. Moreover, it was also found 
that the following three conditions are indispensable for 
the appearance of the nodeless spin-triplet superconduc- 
tivity mediated by phonons: (1) The screening effect is 
sufficiently weak; (2) The short-range Coulomb interac- 
tion is strong; (3) The distance between sheets of the 
open Fermi surface is sufficiently large. Here, conditions 
(1) and (2) are satisfied in (TMTSF)2X, in which it is 
known that the electron states are well-described by the 
tight-binding model, and the layer spacing is large. The 
condition (3) is satisfied by the quasi-one-dimensionality 
(os/6s = 0.468) and the hole density n = 1/4. 
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Fig. 8. Momentum dependences of the order parameters 
^^i^y) = \/ pi.ky) ^(.^Vxiky), ky) of the s, py and d-wave states. 
In these states, A{kx,ky) = A{—kx,ky). 
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Fig. 9. Momentum dependences of the order parameters 
^'^{ky) = \/ p(ky) A(kFx{ky), ky) of thepx and dxy-wave states. 
In these states, A{kx, ky) = —A{—kx, ky). 



In (TMTSF)2X, the pressure-temperature phase dia- 
grams of superconductivity and SDW were obtained. "'^ 
The pressure dependence may be reproduced by taking 
into account the pairing interactions mediated by antifer- 
romagnetic fluctuations, since those interactions have at- 
tractive triplet components.^' We have shown that the 
Px-wave interaction mediated by phonons also decreases 
to some extent as the pressure increases. The explanation 
of the experimental phase diagrams in the present theory 
including the effect of the spin fluctuations remains for 
a future study. 
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